Detecting non-magnetic excitations in quantum magnets 
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Many unconventional quantum phases host special non-magnetic excitations such as photons and visons. We 
discuss two possible ways to detect these excitations experimentally. Firstly, spin-lattice coupling mixes the 
excitations with phonons. The phonon spectral function acquires new features that can be detected by neu- 
tron or X-ray scattering. Secondly, valence-bond fluctuations translate into charge density fluctuations on non- 
bipartite lattices. Such charge fluctuations can be characterized by conventional spectroscopies such as Terahertz 
spectroscopy. Observation of exotic singlet excitations would provide positive identification of unconventional 
quantum phases in frustrated antiferromagnets. 



The search for exotic quantum phases [1] in frustrated anti- 
ferromagnets has been one of the main challenges in the field 
of strongly correlated systems. Such phases are believed to 
emerge when long-range order is destroyed by competing in- 
teractions and strong quantum fluctuations. In a typical situ- 
ation, SU(2) symmetry remains intact down to zero tempera- 
ture so that excitations can be classified by their spin. While 
magnetic excitations can be studied by techniques such as 
inelastic neutron scattering (i.e spinons [2, 3]), exotic non- 
magnetic excitations remain elusive. 

Indeed, novel non-magnetic excitations are predicted in 
generic unconventional magnetic phases. Consider the sem- 
inal resonating valence bond (RVB) phase [4]. Extensive in- 
vestigations of quantum dimer model (QDM) [1] revealed that 
there are two types of dimer liquid in two and three dimen- 
sions. The U(l) liquid exists on bipartite lattices in three 
dimensions. Its low energy excitations are transverse gap- 
less fluctuations of dimer density, or "photons" of an emer- 
gent U(l) gauge theory [5]. On the other hand, the Z2 liq- 
uid appears in non-bipartite lattices and possesses topological 
order The low energy excitations are Z2 vortices, "visons" 
[6, 7]. While a single vison is a non-local object, excitations 
of even number of visons correspond to dimer density fluctu- 
ations [8, 9]. 

While these results were obtained in the QDM, generic con- 
structions [10, 11] exist in which QDM's are low energy limits 
of SU(2) invariant spin models. Such phases and excitations 
could exist in low energy limits of Heisenberg model thanks to 
universality. The authors of an extensive DMRG [12] study of 
spin- 1/2 Heisenberg antiferromagnetic model on the kagome 
lattice concluded that its ground state is a Z2 spin liquid. Stud- 
ies of the multi-spin exchange model [13, 14] on triangular 
lattice found a gapped spin-liquid phase that looks like a Z2 
liquid. Both models are reaUzed in real materials (For reviews, 
see [1] and [15]). Observing singlet excitations in these mate- 
rials would be positive evidence of the existence of Z2 liquid 
phase in nature. 

In this letter, we discuss two general ways to experimen- 
tally probe singlet excitations in quantum antiferromagnets. 
Firstly, singlet excitations mix with optical phonons through 
spin-lattice coupling. For suitable parameters, this leads to 
new features in the phonon spectral function which can be de- 
tected by neutron or X-ray scattering. Spin-lattice coupling 
has been an exciting topic throughout the years. Its study 



was pioneered by the discovery and characterization of spin- 
Peierls transition [16-19]. It was realized [20, 21] that sim- 
ilar mechanisms can induce long range order in the highly 
frustrated Heisenberg antiferromagnet on the pyrochlore lat- 
tice. Magnetoelastic splitting of degenerate optical phonons 
was observed in ZnCr204 [22] and a number of other com- 
pounds. The strongest effect ( 10% splitting) has been seen in 
MnO [23]. Wang and Vishwanath generalized the idea to lo- 
cal phonon [24]. Dynamical effects of phonons were explored 
[25] motivated by spin-Peirels compound CuGeOs [26, 27]. 
These studies demonstrated the importance of both static and 
dynamical effects of the spin-lattice coupling. 

The second way exploits the ability of singlet excitations 
to couple directly to an electric field. Bulaevskii et al [28] 
first discovered that some magnetic ground states and exci- 
tations of certain Mott insulators have nonzero local electric 
charge or current. In particular, fluctuations of valence-bond 
densities induces electric dipoles. The coupling is stronger 
for weak Mott insulators due to smaller U/t where t is the 
hopping amplitude of electrons and U is the onsite repulsion. 
Spin-lattice coupling can lead to the same effect [28]. Con- 
sequently, valence-bond density fluctuations couple to elec- 
tromagnetic radiations directly. Conventional spectroscopic 
techniques may be used to directly detect non-magnetic exci- 
tations. 

The rest of the paper is organized as follows. First, we dis- 
cuss the mixing between optical phonons and singlet excita- 
tions. After introducing a general formulation, we explore its 
consequences for both U(l) and Z2 liquids. We then discuss 
possible charge signature of singlet excitations in the context 
of spin-1/2 Heisenberg antiferromagnetic model on kagome. 
Finally, we conclude our paper by discussing possible discov- 
eries of singlet excitations in real materials. 

Optical phonon: Let us motivate the first mechanism in the 
simplest context. Consider spins interacting via Heisenberg 
exchange whose strength depends on the distance between the 
spins [20-22]: 



J{R + u)S, ■ Sj w J(i?)S, • Sj 



dJ_. 
dr 



\r=R 



{S,-S,)U. (1) 



R is the equilibrium distance between the two spins and u is 
the elongation of the bond, u couples linearly with bond oper- 
ator Si • Sj which measures the singlet density on bond {ij). 
The coupling mixes the singlet excitations and phonons. The 
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phonon spectral function will acquire features of the singlet 
excitations. 

To elaborate on this idea, we consider following Hamilto- 
nian on a general lattice: 



1 



w^u(r)^ + /u(r) • V(r) + 



(2) 

Hs is the spin hamiltonian and u(r) is the displacement of ion 
at site r. We adopt the Einstein phonon model with mass of 
the ion assumed to be 1 for simplicity. Based on the model in 
equation 1, V(r) field is defined as follows [24]: 

"^W^T J2 err'(err'-VrJ(r-r'))Sr-Sr'. (3) 
•' r'GW 

{ r } is the set of neighbors of site r and e^r' = (r — r')/|r — r'|. 

/ = J2r'e{r} ^rr' ' Vr Jr-r' IS the Spiu-latticC COUpliug. 

For small ion displacements, the full phonon Green's func- 
tion in the random phase approximation (RPA) is: 

(r 1 , 1 1 ; ra , ^2 ) = G^^/^ (r i , i i ; ra , ^2 ) -/^Xa^J (i"i , ^ i ; r2 , ^2 ) 

(4) 

where G is the full phonon Green's function, G is the bare one 
and X is the time-ordered bond-bond correlation function: 

1-2,^2) = (T{Fa(ri,ti)y^(r2,i2)}) (5) 
In the Fourier space, relation 4 is written as: 

G-'pik,uj) - G~'pik,uj) - fxap{k,io). (6) 

We explore the consequences of eqn 4 and 6 in two exam- 
ples in QDM's. QDM's [1] are low-energy effective models 
describing the dynamics of nearest neighbor singlets. They 
provide generic playgrounds independent of the underlying 
spin models. 

The first example is U(l) RVB liquid phase on cubic lat- 
tice [5, 29]. The QDM on cubic lattice has two phases [5]: 
the staggered valence bond crystal phase and the U(l) liquid 
phase. The low energy physics of the liquid phase is described 
by the following Hamiltonian in the continuum limit [5]: 



1, 



1 



paB^ +p4(V X B)^ 



(7) 



In the Coulomb gauge = 0, V • A = 0, the electric and 
magnetic fields are expressed as E = 9iA and B = V x A. 
On the lattice, the magnetic field B is defined on the bonds 
[29]: 



-Bo(r) 



na{r) 



(8) 



where Q = (tt, tt, tt) and na(r) is the number of dimers on 
the bond connecting r and r + a {a = x,y,z). z = 6 is the 
coordination number of cubic lattice. 

To include the spin-lattice coupling, we write the bond op- 
erator S(r) • S(r + a) in terms of the gauge field. The bond 



FIG. 1: The new phonon mode around Q is illustrated on a two 
dimensional slice of cubic lattice. The mode is induced by the "pho- 
ton", transverse fluctuation of dimer densities. Bond thicknesses re- 
flect dimer densities. 



operator amounts to two operations on a general dimer cover- 
ing. Its diagonal part counts the number of dimer on the bond 
while the off-diagonal part flips the dimers around plaquette 
to which the bond belongs. The plaquette-flipping operator 
translates into E^ [5]. As a result, the off-diagonal term is ir- 
relevant in the renormalization group sense. The spin-lattice 
Hamiltonian translates into the following compact form in the 
continuum limit: 



H 



sp 



dV/iB • u 



(9) 



where u = e'*^ ''u and /i is the spin-lattice coupling. The 
Hamiltonian for phonon is: 



1 



(10) 



The total Hamiltonian is H = Hp + Hgp + Hg. 

We focus on the phonon spectrum in U(l) liquid phase 
where p4 term can be neglected. Applying equations 4 and 6, 
phonon develops two transverse sound modes around momen- 
tum Q. These modes manifest themselves as new low-energy 
poles in the phonon spectral function. For momentum Q + k 
(k ^ 1), the energy is approximately: 



I P2 



(11) 



The spectrum weight of the modes is proportional to 
fik'^/{ujQ — p2k^)'^. In contrast, the spectrum of the longi- 
tudinal phonon remains unchanged. This is a reflection of the 
transverse nature of gauge fluctuations. Such sound modes 
generally exist in QDM's on other three-dimensional bipartite 
lattices. 

If the spin-lattice coupling is so large that ff > p2CdQ, the 
lattice distorts and the magnetic fluxes condenses. This is the 
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FIG. 2: One-vison and two-vison excitations. A vison at i is created 
by operators on all thie bonds to the left of i. A two-vison excitation 
is related to the dimer density on the bond I between them. 

analog of spin-Peierls transition. We will discuss it in a future 
publication. 

The second generic dimer liquid phase is the Z2 liquid 
phase on two and three dimensional non-bipartite lattices [ 1 ] . 
Such a state preserves all lattice symmetries and has a gap to 
all excitations. The system possesses topological order Con- 
sider a Z2 liquid on a cyUnder, the state belongs to the even 
or odd topological sectors if a cut around the cylinder crosses 
even or odd number of dimers. The low energy excitations 
are visons [6, 7], Z2 vortices residing on the sites of the dual 
lattices. At the Rokhsar-Kivelson (RK) [30] point, the ground 
state is an equal amplitude combination of all possible dimer 
state. To obtain a vison at site r, we make a cut [ 1 ] from the 
site to the lattice boundary (Fig 2). A state with one vison is 
a combination of all possible dimer coverings. However, the 
amplitude of a dimer covering is negative if the cut crosses an 
odd number of dimers. Visons are nonlocal objects. A local 
operator such as Sr • ^^+1 involves even number of visons. 

We consider two cases. In the first scenario, two-vison ex- 
citations have definite energy and momentum. This makes 
sense if visons form bound pairs. Such a bound state has been 
observed [31] in numerical studies of QDM on triangular lat- 
tice. It is also possible that the vison dispersion is almost flat 
[32-34] so that the two-vison continuum is very narrow. In 
the former case, the two-vison excitation induces a new pole 
in the phonon spectrum in the same fashion as the U(l) case. 
For the latter, a very narrow continuum will be present in the 
low-energy spectral function of phonons. The spectral weight 
of either the pole or the narrow continuum is determined by 

/u!q where uj^ is the energy of the two-vison excitation. 

If the two-vison excitations form a broad continuum, the vi- 
son has a large bandwidth with N minima in the first Brillouin 
zone. We introduce coarse-grained field operators around 
these minima: 

Mr) = 4^Y. ^^"".qC-ic'^J- (12) 

where a — 1 . . .N and A is the momentum cut-off set by the 
lattice scale. In terms of the field operators, the vison-phonon 
coupUng can be written in the following form: 

+/Se"*^*'°"*'^V;(r)</>^(r) + h.c) . (13) 




FIG. 3: A dimer-covering state of the kagome lattice. A vacuum 
triangle carries an electric dipole illustrated by the black arrow as- 
suming Q > 0. Vertices around the defect triangle carry charge — Q 
while other vertices carry Q. 



g jp and /^^ can be derived from a microscopic model. They 
can also be fixed up to an overall constant by s ymmetries. The 
energy of vison fields is cjq = -^/m^ + v'^q^. For simplicity, 
we assume jtIq, = m and = v. Since visons are gapped, 
only 4)a4>p and its complex conjugate contribute to x- 

As u) approaches the two-vison continuum from below, 
x(w,kQ + k^) diverges logarithmically at w = 2m. This 
leads to appearance of a new pole in phonon Green's function 
4. The residue of the new mode increases as the pole moves 
away from 2m. For appropriate parameters, such modes can 
be observed. 

As 771 decreases, the new phonon modes are pushed toward 
zero energy. At a critical m, one (or several) of them con- 
denses. The system becomes a valence bond crystal. This 
is the analog of Spin-Peierls transition in Z2 liquid. We will 
discuss it in detail in a future work. 

Charge signature: On non-bipartite lattices, singlet density 
fluctuations generate electric dipoles [28]. We consider three 
spins on a equal lateral triangle interacting antiferromagneti- 
cally. The exchange energy is minimized by combining two 
of them into a singlet. The induced electric dipole lies in the 
plane of the triangle normally to the singlet bond [28] (Fig- 
ure 3). In other words, charge 2Q accumulates on the free 
spin while the two spins forming the singlet carry charge —Q 
each. If the dipole is induced by higher-order perturbations in 
a weak Mott insulator, Q > is proportional to (t/U)^. The 
sign and the magnitude of spin-lattice coupling determines Q 
when the dipole is generated by magnetostriction. 

Consider spin- 1/2 Heisenberg antiferromagnetic model on 
the kagome lattice (Figure 3), a network of corner shar- 
ing triangles. The low energy states are dimer-covering 
states with maximum number of nearest neighbor singlets. 
A quarter of triangles, the so-called "defect triangles", lack 
singlets[35, 36]. The uneven distribution of singlet densi- 
tiestranslates into inhomogeneous charge densities. A simple 
counting shows that the vertices of defect triangles each carry 
—Q while all other vertices have Q. 

An applied AC electric field can be used to induce the mo- 
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tion of defect triangles. For typical sample size (L = 1 mm) 
and singlet excitation energy (cj = 0.1 meV), the applied field 
is approximately uniform, q = ui/c ^ 1/L. The scattering 
of the applied field will provide information about the singlet 
spectrum of the system at q = 0. Suitable techniques include, 
for example. Terahertz spectroscopy. 

Discussion: In this letter we described how exotic singlet 
excitations can be detected, at least in principle, by exist- 
ing spectroscopic methods. While focusing on specific mod- 
els, we stress that the two mechanisms described are model- 
independent as long as the global SU(2) symmetry is intact. 

Currently two classes of materials could host the Z2 liq- 
uid state. The first set of materials including ZnCu3(OH)6Cl2 
(herbertsmithite), Cu3V207(OH)2 -21120 (volborithite) and 
BaCu3V208(OH)2 (vesignieite) (See [1] and [15] for reviews) 
realize the 5 = 1/2 Heisenberg antiferromagnetic model on 
the kagome lattice. Yan et al [12] presented some evidence 
that the ground state of the model is a Z2 spin liquid. The vi- 
son spectrum were studied by several groups [32-34]. They 
identified the location of low energy singlet excitations in the 
reciprocal space. Recently, single-crystal sample of herbert- 
smithite was synthesized [37]. Measuring the spectral func- 
tion of phonons in the magnetic energy range (up to a few 
meV) could reveal novel singlet excitations. It would also be 
interesting to measure its spectrum at q = using conven- 
tional spectroscopies. 

The second class of materials including k-(BEDT- 



TTF)2Cu2(CN)3 and EtMe3Sb[Pd(dmit)2]2 (See [15] for a 
review) realizes the multi-spin exchange model on triangu- 
lar lattice. Studies [13, 14] show that the ground state of the 
model can be a gapless or a gapped spin liquid for different 
parameters. The gapped liquid phase resembles the Z2 liquid 
phase with a large number of singlet excitations within the 
spin gap [13]. While k-(BEDT-TTF)2Cu2(CN)3 [38] seems 
to host the gapless liquid phase, the flexibility of the material 
family k-(BEDT-TTF)2X raises the hope that the gapped liq- 
uid phase is the ground state for some other member whose 
singlet excitations can be observed by studying the phonon 
spectrum. 

Beyond quantum magnetism, singlet excitations are be- 
lieved to be important for other strongly correlated systems 
such as high temperature superconductors [39]. We speculate 
similar couplings between optical phonons and singlet exci- 
tations also exist. It would be very interesting to search for 
the trace of singlet excitations in phonon spectrum in these 
systems. 
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